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Abstract 

An extended Hubbard model with phonons is considered, q-coherent states relative to 
the superconducting quantum symmetry of the model are constructed and their properties 
studied. It is shown that they can have energy expectation lower than eigenstates constructed 
via conventional processes and that they exhibit ODLRO. 
PACS: 74.20.-z, 05.30.Fk, 71.27+a 

The Hubbard model jTJ, despite its simplicity, is certainly one of the important tools in 
condensed matter physics. It is one of the simplest models to describe itinerant interacting 
electrons on a lattice (for review see e.g. It is also believed that it may be used for the 

description of high T c -superconductivity. 

In [Sj, the so-called 77-pairing mechanism of superconductivity was introduced. This allows 
the construction of states exhibiting off-diagonal long range order (ODLRO) which is regarded 
as a definition for superconductivity |l].The n operators involved in this construction generate 
an su(2) algebra. More precisely, it has been shown ^ that the standard Hubbard model 
have a SO (4) = SU(2) x SU(2)/Z2 symmetry. The first SU(2) symmetry (spin symmetry) 
is useful for the description of the antiferromagnetic properties of the electron systems. The 
second SU(2) symmetry (pseudo spin), when broken, is responsible of the superconductivity 
of the model. 

In an extended Hubbard model with phonons was shown to violate the superconduct- 
ing SU(2) symmetry and that it can be only restored as a quantum (deformed) [7 9 (sit(2)) 
symmetry. 

In this letter, motivated by an idea [7J that coherent states |81l51llU| may be very relevant 
for the study of the standard Hubbard model, we introduce q-coherent states related to the 
extended Hubbard model of These states are nothing but the coherent states related to 
the quantum algebra [/ 9 (sm(2)) (see e.g. |11|1 and are shown to have, under some conditions, 
energy expectation lower than states constructed through the ^-pairing mechanism and even 
from the Hamiltonian eigenstates constructed in 0. We also show that this q-coherent states 
exhibit ODLRO and are thus superconducting. 

First, we begin by some preliminaries about the Hubbard model and recall the extended 
one. 
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Consider a D-dimensional lattice, A, with L sites. The standard Hubbard Hamiltonian 
pQ describing a system of itinerant interacting electrons in A, written in a grand-canonical 
formalism, is 



where 



( - n(rii^ +TH,i) + u rn,xni,ij, (2) 

Hir n ~ loc) = \t J2 E (3) 



ieA 



2 

<i,j> ct=T4 

We have used the conventional notations, namely, 

• ai,a and a^" CT are the canonical Fermi operators describing electrons on the lattice A and 
obey the usual anticommutation relations. 

• i and j are (real space) sites in A and < i,j > means that the summation is carried over 
nearest neighbor sites. 

• a =1, | denoting spin down and spin up respectively 

• Hj.CT = afa-Cii^cr denotes the number operator of electrons with spin a at site i. 

• fi is the chemical potential, u the on site repulsion energy and t the hopping amplitude. 

The model described by (0 has a natural 5(7(2) x SU(2)/Z2 symmetry ]5 at half filling. 
The first SU(2) symmetry, usually called spin or magnetic symmetry, is generated by the 
following spin operators 

S + =E<T ffl *-l> S- = (S + y, S* = 5^(71^-71^) (4) 

ieA ieA 

and is reminiscent of the antiferromagnetic properties of the system, see e.g. 0. 
The second SU(2) symmetry is generated by the 77 — pairing operators 

V + =J2<l<1 ' r > = (l + ?> rf = ±(N-L) t (5) 

ieA 

where N = S nt := S ' S ' 7ii, CT = ^ ' riij + rii,j, denotes the total number operator. 
ieA ieA o=T4 ieA 

It is this latter symmetry that, when broken, leads to superconducting repaired eigenstates 
of the Hamiltonian (for this reason it is called a superconducting symmetry). 

Both symmetries can still be conserved when some additional terms and interactions are 
added to the Hamiltonian see e.g. |12l I13I FRl I15I I16| . However this is not the case in 
It is indeed, the latter (superconducting) SU(2) symmetry that is modified (deformed) there. 
Their reasoning is the following. 

Montorsi and Rasetti considered the extension of the Hubbard model, described by Q, 
by including phonons in it to describe the ion vibrations. The relevant Hubbard Hamiltonian 
then is 

H = H ( e l ° c) + H { p l ° c) + H$$ h + H (non - loc) . (6) 
H^°°^ is given in the other terms are 

l P h 2M 



H ~ l ° c) = E m + \ Muj2x - (7) 



H^ ph = -A^Kt+^i)^, (8) 

ieA 

H(non - loc) = (U^a^ + h.c), (9) 

<i,j> ff=T4 
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where the hopping amplitude depends now on the ion displacement and is given by 

Uj = t exp |c(£j — Zj) j exp |ifc(pi — Pj)|. (10) 

Pi and a;; are local momentum and displacement operators for an ion on site i, [xi,pj] = ihSij, 
and are both assumed to commute with the fermionic operators. Q gives the kinetic term 
for phonons, considered as a set of decoupled Einstein oscillators with frequency ui and mass 
M. Finally, A is the coupling constant for the local phonon-electron interaction (|HJ. 

The important fact to notice here is that when A / the magnetic symmetry (|1J is 
conserved while the superconducting one is broken. It can be restored only as a quantum 
(deformed) symmetry, i.e. [/ 9 (sit(2)) , as proven in [S) (see also |17|1. the reasoning being 

basically the same as the one in [Sj- Namely, first a local U q (sm(2)) symmetry (i.e. a symmetry 
for the local part of the Hamiltonian) is constructed, it is generated by the following operators 

^+)=e-^a+T<|, Kt^iK^y, = i(r»,, T + «,,;- l) . (11) 

For the moment <f> is an arbitrary real parameter. The operators (111! obey the usual U q (sit(2)) 
commutation relations (see e.g. |18l 

[K^,K^] = ±K^ , [k\ + \k\-\ = [2Kl% , (12) 
with q — e _Q! , a being an arbitrary complex. The box function \] q is defined by 

[X] q = q — V - • (13) 

q — q 1 smh(a) 

For latter use we also define the q-factorial function as [n] 9 ! = [n] 9 [n— l] q ...[l] 9 , [0] 9 ! = 1 for 
n being an integer. 

The requirement that the operators Q110 commute with the local Hamiltonian J/^° c ' yields 
the following constraints on the parameters involved 

If 2A \ 2A 
»=2\ U -M^)> ^ = KM^- (14) 

The first constraint being the analogue of the half-filling requirement in the standard Hubbard 
model (i.e. A = 0). The second one constraining the choice of the parameter (f>. 

The extension of this local symmetry to a global one over A is achieved by means of the 
corresponding coproduct 0. Note that, to be able to do this, it is necessary to adopt some 
ordering of lattice sites. 

The resultant (global) U q (su(2)} algebra is generated by 

K { z) = J2 K i Z) ( 15 ) 

K (+) = ^-yoj jj e -a-*<*> (16) 

j£A k<j k>j 

K (-) = ( Ji -(+))t < (17) 

The phase factor e lG 'j, with G = (ir, ir), will be useful in what follows. 

The above quantum algebra [7q(sw(2)) describes a global symmetry of the model (i.e. the 
generators I15ll7t commute with the total Hamiltonian JfjJ if and only if 

Re(a) = ^ and k = 2cj>, (18) 

in addition to these two constraints, in |17| it has been shown that there are also additional 
constraints on the ordering of the lattice sites for the symmetry to hold. The authors of |17| 
thus concluded that this symmetry holds only on one-dimensional lattices. 
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The first relation in 11811 gives a constraint on the real part of a, so that its complex part 
can be chosen arbitrarily, without loss of generality we will take it null. The second relation, 
mainly, means that we can interpret ^ as the parameter of a Lang-Firsov transformation |17|. 

In their paper 6 Montorsi and Rasetti also constructed eigenstates for the Hamiltonian 
© with energy expectation lower (or at least equal) to that of states constructed via the 
conventional ^-pairing mechanism. These states are constructed from the vacuum state by 
successive action of the operator _K^ + ': 

\4> n > = f3(n, L) (K w ) n \vac > , (19) 

where the normalization constant is given by 

«-«-(&&)*■ * 

This is the case when the t/ ? (.su(2)) symmetry holds. 

When the model parameters have values such that this symmetry does not hold any more, 
one can still prove the following relation 

[H,K {±) ] = ±EK (±) , (21) 

where 

E = u -^~ A w^- (22) 



Now that we have reminded all the necessar materials, we can start developing our idea. 
This is an adaptation of an idea anticipated by Solomon and Penson in for a standard 
Hubbard model (SU(2) symmetry). In fact, they have used coherent (pairing) states in study- 
ing the Hubbard model and proved that this states can be useful for the description of the 
model. In the following, we shall follow a similar reasoning and use q-coherent states for the 
description of Montorsi and Rasetti's Hubbard model. 

Since this model have [7 9 (sit(2)) symmetry it is natural to use q-coherent states, for in- 
stance see |11| and references therein. 

The definition we adopt for this states is the following 

\v>= N~^(\v\ 2 ) exp^vK^) \4> L > . (23) 
Note that we have used the maximal state 

\4>l>= j^^hO, (24) 
and the deformed exponential function 



The normalization constant in 1231 can easily be computed using some standard q-series 
technics |20) and is given by 

AA(M 2 ) = (l + W\ 2 ) L q (26) 
:= f[(l + o L - 2k+1 H 2 ) 

I |2n 
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where the q-binomial is defined by 



(27) 



The states \v > are clearly not eigenstates of the Hamiltonian JSJ since they involve 
eigenstates (\<j> n >) with different n's. They posses, however, many interesting properties 
that we shall derive in the following. 



First of all, let us notice that in the limit L — > oo, \u > is an eigenstate of 



Thus, 



it (\i> >) can be seen as a q-harmonic oscillator coherent state, see e.g. |21|. This is not 
surprising since in the standard Hubbard model the coherent states constructed in '7J behave 
similarly, i.e. in that limit their coherent state is seen as a (standard) harmonic oscillator 
coherent state. 

Let us evaluate the energy expectation in a q-coherent state 123H . For this purpose we will 
use the following formulae, derived from 12 1 1 - 1221 . 



[H,(K (±) ) n ] = ±nE(K (±) ) n , 
[H,exp {vK<~- ) )] = -vEK^expJvK^) 



We obtain the following formula for the desired expectation value 
<v\H\v>=M-\\v\ 2 ) E V|i/| 2 " (L 

L n J q 

n=0 

or, equivalently, 



< v\H\v >= LE -Af^diy] 2 ) E 



n=0 



(28) 



(29) 



(30) 



In contrast with the results of [7], this formula for the energy expectation can not be 
further simplified. However, it is easy to see its behavior as the different parameters involved 
change. For instance, in the limit q — > 1, as expected, one recovers a similar formulae to that 
obtained in |7j. 

Also, using the fact that [n] q > n is always true, it is easy to see that as q tends to 1 the 
energy expectation 1291 tends to E and after a critical value of q it gets even smaller than it. 

The same reasoning on \v\ implies that the energy expectation decreases as \v\ gets bigger 
and it gets smaller than E after a critical value of \v\. 

All this is confirmed by the plots (see plot a, plot b and plot c) 



Another important feature of the q-coherent states 1231 1 is that they exhibit off-diagonal 
long range order (ODLRO) pi] and thus are superconducting states. In fact, the relevant 
off-diagonal matrix element of the reduced density matrix p2, by considering the states \<j} n > 
11911 . is equal to < (p n \cij^aj , a r ,ia r ,-\\4>n >■ And it was found |S] to be equal to 



e iG(r-s) e a(JV+l-|r-s|) 



L 
n—1 



r L 



7^0 



(31) 



Using this result and I12MI . one can show that the relevant expectation value to be evaluated, 
< u\a^ ,a r ,io,r,t \ v >j equals 



e iG(r-s) e a(N+l-\r-s\) f 1 + M 



/f-H\v\>) 



(32) 
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which is also different from zero (for large values of \r — s|). Thus the q-coherent states \u > 
are superconducting. 

As a matter of fact, this is not surprising . It has been argued in |13l 1151 that supercon- 
ductivity based on ^-pairing is a generic rather than an exotic phenomenon. On the other 
hand the states \<f> n > 11911 were shown to have pairing Moreover, the coherent states 
of 7 , which were constructed using an 77-paired state, posses ODLRO. All these, permit us 
to conclude that superconductivity of coherent states, constructed from states which possess 
pairing, is not an exotic phenomenon. 

In summary, we have constructed q-coherent states for the superconducting [7 9 (sm(2)) 
symmetry of an extended Hubbard model with phonons. We have shown that they have energy 
expectation lower than that of the eigenstates \cj> n > constructed in It was also shown that 
these states exhibit ODLRO, thus are superconducting. We concluded that these property 
(ODLRO) should be expected for any coherent states constructed using pairing states. 
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Plot c. 

Energy expectation < u\H\u >, eq(29), versus: 

- \v\ for q — 0, 9, L — 10 and normalized E; Plot a. 

- q for \i>\ =3, L — 10 and normalized E; Plot b. 

- |^| and q for L = 10 and normalized E; Plotc. 
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